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Abstract 

A simple model is constructed which allows to compute modified dispersion relations with effects 
from loop quantum gravity. Different quantization choices can be realized and their effects on the 
order of corrections studied explicitly. A comparison with more involved semiclassical techniques 
shows that there is agreement even at a quantitative level. 

Furthermore, by contrasting Hamiltonian and Lagrangian descriptions we show that possible 
Lorentz symmetry violations may be blurred as an artifact of the approximation scheme. Whether 
this is the case in a purely Hamiltonian analysis can be resolved by an improvement in the effective 
semiclassical analysis. 
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I. INTRODUCTION 



One of the cases where observations of quantum gravity effects have been imagined relies 
on modified dispersion relations for matter (such as photons or neutrinos) travelling on a 
quantum gravitational background [lj. Planck scale effects are expected to be negligible 
in standard circumstances, but observations of highly energetic particles travelling long 
distances or other high precision experiments may set bounds on possible effects. Such 
effects have been interpreted as indicating a violation of (standard) Lorentz symmetry 
for which, as expected, the observational bounds are rather stringent Q]. The theoretical 
foundation, on the other hand, is still open and actively debated. For instance, while trying 
to link the breaking of Lorentz symmetry to a privileged frame has been argued to be in 
conflict with our current understanding of field theory one can try to invoke a deformed, 
rather than broken, symmetry in the form of Doubly Special Relativity |5( . Independently 
of these field theoretical considerations, the task of candidate quantum theories of gravity 
is to provide reliable estimates on the magnitude of expected modifications to the standard 
dispersion relations to be compared with observations. 

n 

One such candidate is loop quantum gravity p which leads to a discrete structure of the 
geometry of space. This discreteness can be expected to lead to small-scale corrections of 
dispersion relations, just as the atomic structure of matter modifies continuum dispersion 
relations once the wave length becomes comparable to the lattice size. There have been 
several studies already which derive modified dispersion relations motivated from particular 
properties of loop quantum gravity , but at this stage the control on the calculations 

is insufficient. The difficulty lies in the fact that loop quantum gravity is very successful in 
providing a completely non-perturbative and background independent quantization of gen- 
eral relativity which makes it harder to re-introduce a background such as Minkowski space 
over which a perturbation expansion could be performed. Techniques for constructing semi- 
classical states are available and still being developed further [ly, [l]] , but the calculations 
toward modified dispersion relations are very complicated. Moreover, the answer cannot be 
expected to be unique but to depend on several parameters as well as quantization choices 
in the full quantum theory. 

In the first part of this paper we develop and study a simple model which allows us to 
introduce crucial properties of loop quantum gravity into the Hamiltonian of a matter field. 
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As we will show, the model captures essentially all the effects that have been considered 
so far in loop motivated calculations of modified dispersion relations even at a quantitative 
level. We can also see how different quantizations would change the results, and which 
quantization choices should have the largest effect on the order at which corrections occur. 
Thus, we have the freedom to change basic objects according to the possibilities of loop 
quantum gravity, but a much simpler and more immediate way to check the consequences. 

The essential idea in constructing the model is to consider space as being made of ho- 
mogeneous patches defining a lattice on which the matter Hamiltonian, in particular its 
space derivatives, will be discretized. This models the discrete structure of loop quantum 
gravity, but could also be used classically as an approximation (the metric field is then sim- 
ply considered as a piecewise constant rather than continuous function). Such a classical 
approximation would become better and better if we choose smaller and smaller patch sizes. 
A second ingredient from quantum geometry then is that the patch geometry must be quan- 
tized (which is readily done for homogeneous or even isotropic patches 12j)- This implies 
additional, quantum geometric corrections which grow with shrinking patches. Thus, with 
effects from quantum geometry there is a non-zero patch size leading to a minimal deviation 
of the effective matter Hamiltonian from the classical one. 

This model can be formulated at different levels of complexity which allows to consider 
more realistic situations and also to bring it closer to what one would have in full loop 
quantum gravity. In this paper, we only consider the simplest construction using isotropic 
patches of equal size, and only couple a scalar matter field with simple discretizations. Still, 
as we will show below, this simple model captures essentially all the effects that have been 
considered so far in loop motivated calculations of modified dispersion relations even at a 
quantitative level. Although there are still gaps between the model and full loop quantum 
gravity, it gives - because of its simplicity - a more direct link between quantum gravity 
phenomenology and the full theory. 

In the second part of the paper, we study the possibility that the apparent Lorentz viola- 
tion of dispersion relations obtained in our simple model, as well as in the more sophisticated 
treatments e a result of the approximation scheme rather than of the theory itself. 

Common to all these computations is that they derive corrections to the matter Hamiltonian. 
This is natural in the setting of loop quantum gravity because the quantization of spatial 
geometry is readily available whereas 4d covariant quantities are harder to quantize. On 
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the other hand, by way of examples we will demonstrate that a purely Hamiltonian analysis 
is much more subtle than a Lagrangian one, and discuss how a perturbative Hamiltonian 
analysis can be improved in order to draw reliable conclusions. 

At first it might seem that Hamiltonian and Lagrangian descriptions are completely 
equivalent, which certainly is the case when theories are considered exactly. However, as we 
will demonstrate, the situation changes when approximation schemes, in particular pertur- 
bation theory, are employed. When higher derivative corrections are involved, the Legendre 
transform does not commute with expansion in a perturbation series such that the perturbed 
Hamiltonian and Lagrangian formulations are not necessarily equivalent. This is in par- 
ticular the case for theories non-local in time, in particular when time is not a continuous 
parameter but 'discrete'. This might well be the case in quantum gravity and specifically 
in loop quantum gravity . We will show in examples that going over to an approximate 
continuum Lagrangian and then to the Hamiltonian description will lead to a Hamiltonian 
that one could not have obtained from a Hamiltonian model with continuous time in the 
spirit of what has been proposed for loop quantum gravity. Our conclusion will therefore be 
that many of the calculations done up to now (including the first part of this paper!) can 
only yield preliminary results and that a definite answer to the question of Lorentz violation 
by loop quantum gravity will have to await a more complete treatment. We will conclude 
with an outlook on possible strategies in this direction. 



II. QUANTUM CORRECTIONS TO THE SCALAR FIELD HAMILTONIAN 

The simplest way to couple loop quantum gravity (LQG) to a free scalar field is via its 
Hamiltonian 

H = Jj 3 x QiV((det q)- l/2 pl + Vd^q ab d a <H0)) (1) 

where £ is a Cauchy surface of the space-time manifold M (assumed to be globally hyper- 
bolic). The more complete treatments in the literature proceed by quantizing the gravita- 



tional part of this Hamiltonian with LQG methods 3| and then take expectation values in 
a semiclassical state. That state comes with a discretization of the spacial slice £ and as 
a consequence, the partial derivatives in the classical expression Q are changed to lattice 
derivatives. Further differences to the classical expression result from quantum corrections 
to the classical values of (det <2 , ) ±1/ ' 2 , q ab upon taking expectation values of the corresponding 
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operators in the semiclassical state. 

To compare the corrected Hamiltonian obtained by this method with the standard ex- 
pression (|T|). the lattice derivatives in the former are expanded in a Taylor series to obtain, 
by collecting the lower orders of this expansion and of the other quantum corrections, an 
effective Hamiltonian. Plane waves are solutions to the equations of motion generated by 
this effective Hamiltonian and the corresponding dispersion relations contain corrections 
compared to the standard one. We emphasize that in this section we follow the standard 

nn 

procedure jg, |9( . 

To compute the corrections for the Hamiltonian in a simplified way, let us now propose a 
model which includes the expected properties (and which can always be made more compli- 
cated to be more realistic). Let us discretize space into patches on which we can assume the 
geometry to be approximately isotropic. Each patch a then carries two real numbers, one for 
the scalar <p a and one for its momentum p^ a , and an isotropic semiclassical quantum state 
ip a for the geometry. (There is also a lapse function, a real number, per patch which is not 
so important for our purposes.) This corresponds to a scalar on a classical geometry made of 
patches of a size given by the expectation value of the volume operator in the semiclassical 
state. 

Again, corrections in the Hamiltonian are of two kinds: Since it contains space derivatives 
of </>, which have to be replaced by finite differences, there is a discretization error. In 



addition, geometric quantities like vdet q and its inverse have to be replaced by expectation 
values of the corresponding operators. Choosing each patch to have coordinate volume one 
(otherwise, there would be unnecessary coordinate factors), one obtains 

tfdisc = lj2 N °[pl aiidetq) ~ 1/2) ^ 

a ' (2) 
+ ( E i E i/ V / detg)^ a (0a+e / - <p a - ei ){(t> a +ej ~ 0a- ej )/4 . 

Here, a subscript ip a means taking the expectation value in the state if) a , and a + ej denotes 
the neighboring patch in direction /. For isotropic patches, the expression simplifies to 



I 

,2 • 



(3) 



where p = a is the isotropic densitized triad component. 

So far, there are many parameters to specify the background for the scalar: For each 
patch we have a state, which is characterized by its expectation value for p and its spread. 



At first, one can assume that all patches have the same values, which still leaves us with two 
scales in addition to the Planck length and a wave length. Since the difference corrections 
increase with the size of the patches while the corrections for inverse powers increase with 
decreasing size, the first scale, the scale factor a of the isotropic patches, can be fixed by 
requiring a minimal sum of those corrections. 

Specifically, we can relate the discrete scalar field values <p a to a continuous field 4>(x) by 
4> a = 4>(x a ) where x a is a point in the center of patch a. Expanding <p(x) we get 

0a+ej — (fi(x a + e[ ) = <f>{X a ) + (la^e, — X a ) a d a (f)(x a ) 

+ 2^- Xa+e i ~ Xa ^ \ x °t+ei ~ Xa ) 9 a db(p(x a ) + . . . 

Now let us assume that the x a are the vertices of a regular cubic lattice aligned with the 
coordinates on S. Then simplifies to 

<p a+ei = 4>(x a ) + <9/0(x a ) + ^d](p(x a ) + .... (5) 

The squared differences in the Hamiltonian Q are thus approximated by 

~ ^(0 a+e/ - K-etf = + +■■] M , (6) 

/ / 

i.e. we get the second derivative term we need plus higher derivative corrections. Let us 
write 

B :=\Y,di<PdU (7) 
6 i 

for the first order correction. It is certainly not rotation invariant, the symmetry having been 
broken by the introduction of the regular lattice of patches. In a more realistic calculation 
one would work with a random lattice or average over regular lattices with different orien- 
tations to define the semiclassical state. As we are only interested in an order of magnitude 
calculation, we disregard this issue here. 

For the corrections of the momentum term in (j3J) we can use earlier calculations for the 
inverse scale factor resulting in (p _3 ^ 2 ),/> a = o~ 3 + Aa -3 with a quantum correction Aa~ 3 
which for larger a is perturbative in the Planck length. In a triad eigenstate Q], those 
corrections would be p^Aa~ 3 = cp^a~ 3 -£p/a A = ca^N^cf 2 ^ (since = Na 3 (p) with some 
constant c which can be computed once we make a choice on the explicit quantization, while 
a coherent state would result in <i 2 £|/a 4 instead of £p/a 4 . (Thus, the Planck length 
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would be replaced by its geometric mean with the spreading scale d of the coherent state.) 
Using a coherent state thus makes the correction smaller, as expected, but we cannot yet 
tell the order in £p since a is not fixed. To do this, we minimize the total correction to the 
Hamiltonian (J3J) 

ca^N' 2 ^ 2 ^ + aB (8) 



with respect to a. This gives 



for a triad eigenstate and 



x 1/2 

2 tp B 



(9) 



c j2 «2 N- 2 4> 2 



1/2 



a= (10) 
for a coherent state. 

A classical wave solution with our notation (e.g. dimensionless coordinates) has the form 

= exp(i(ak ■ x + Nut)) (11) 

if all patches have the same size a (otherwise, we would have to sum to get the physical 
distance in the argument). Thus, we have an implicit expression for a 

a = ^/(c/2)4cj 2 A 4 /(16vr 4 c'a 4 ) (12) 

where we have expressed B as c'(2iia/ A) 4 2 and d is between 1/9 and 1/3, depending on the 
direction of propagation (this combines the factor 1/3 in B with another factor of \k\~ A Yli kj 
which can be seen to be bound by 1/3|A;| 4 < Y^,i ^ l^| 4 )- Thus we get 

a = ((c/(327rV)4^ 2 A 4 ) 1/6 « ((c/(8vrV))4A 2 ) 1/6 (13) 

where we used ui ~ 27rA _1 which is only approximately true since we expect corrections to 
the dispersion relation (but corrections here would only affect higher order terms). 

Thus, the patch size is a weighted mean of the Planck length and the wave length, with 
a large weight on the Planck length. This means that we need rather small patches, and it 
strongly reduces the order of expected corrections: From the inverse powers of p we expect 
corrections of the order 



4/a 4 oc (£ P /A) 4/3 (triad eigenstate) oc (£ P /Vd\) 8/3 (coh. state) . (14) 
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The order of the corrections from the y/p-term on the other hand is 

a 2 /X 2 oc (£ P /A) 4/3 (triad eigenstate) oc (v^/A) 4/3 (coh. state) . (15) 



These results square very nicely with (t| : First of all in both cases there is a spatial discretiza- 
tion, leading to discretized derivatives and, consequently, higher order derivative corrections 
in the effective Hamiltonian. The characteristic size of the discretization e in j2j is (as a 
here) a weighted mean e m £gA 1_Q . While a was not uniquely fixed there, a = 2/3 would 
reproduce the result here. Also, the order of the first correction due to the higher derivative 
terms was, exactly as in our case, found to be e 2 /A 2 . 

Next, since (of also works with coherent states, there is a parameter (called a there) 
corresponding closely to our parameter d, distributing the width of the state between con- 
figuration and momentum degrees of freedom. There, this parameter is chosen macroscopic. 
One can probably understand this from our result here that the correction (pPfj) decreases 
while the correction (|15jl increases with increasing d. However, the decrease is governed by 
a higher power of d such that larger d are favored. 

Finally, the relative order of the correction due to quantum effects found in joj is 
(e/A) 2 / a_1 . Again, this corresponds precisely to our (£p/a) 4 for a = 2/3. As for comparison 
with the dispersion relations in |s| the analysis proceeds similarly. 



III. THE ISSUE OF LORENTZ IN VARIANCE 



A case in which modified dispersion relations have good chance of being tested is when 
they break Lorentz invariance. This allows correction terms of the form £pE which can 
be high enough for sufficiently large energy E. Possible Lorentz invariant corrections, on 
the other hand, can at most be of the order Ipm with the fixed and limited mass m. 
Accordingly, except for trivially modified dispersion relations that have been discussed for 
quantum gravity phenomenology, all break Lorentz invariance. 

Loop quantum gravity in particular is a Hamiltonian formalism where Lorentz invariance 
is not manifest. (The supposedly covariant twin of loop quantum gravity, spin foams, is under 
much less control currently and not yet suitable for explicit applications; moreover, anomaly 
free formulations may even loose manifest covariance Q].) If we first consider only the 
spatial aspects, rotational invariance is not manifestly broken by the discrete structure since 
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one does not restrict the theory to a fixed lattice. Nevertheless, calculations of dispersion 
relations choose a fixed graph which cannot be rotationally invariant, but do the calculations 
in such a way that the result is rotationally symmetric. So a priori discreteness does not 
necessarily imply violations of symmetries in the approximate classical expressions. For 
Lorentz invariance, however, this is much more difficult to achieve since time does not 
appear directly in the theory. 

The Hamiltonian formulation requires calculations to be based on a lattice in space such 
that only the spatial geometry is manifestly discrete. Nevertheless, dispersion relations are 
computed from classical field equations which involve coordinate time. This time parameter 
is introduced by computing the perturbative matter Hamiltonian on the lattice, and then 
treating it as the Hamiltonian of a classical field theory. Time then appears via the Hamil- 
tonian equations of motion, but only at the classical level. In particular, time is always 
continuous in this setting unlike space, whose discrete structure is responsible for the very 
effects to be computed. 

That the situation in loop quantum gravity is indeed such that the calculations done so 
far introduce Lorentz violations not coming from the theory is suggested by an additional 
complication for this kind of question caused by discrete theories. Discrete theories are non- 
local which implies that they have effective formulations of higher derivative type (when 
differences are to be expanded in a Taylor series to arrive at an effective Hamiltonian or 
action). For higher derivative theories, in turn, the Legendre transform does not commute 
with a perturbation expansion: If we start with a higher derivative Lagrangian in which the 
higher derivative terms can be treated as perturbations, and Legendre transform, then the 
resulting Hamiltonian will not be analytic in the perturbation parameters 1^ . If, on the 
other hand, we first Legendre transform the full Lagrangian then the perturbation expansion 
of the resulting Hamiltonian must obviously be analytic in the perturbation parameter. 

While a Lagrangian formulation would immediately show whether or not Lorentz in- 
variance is broken by correction terms, the Hamiltonian formulation is more indirect. Since 
perturbation and Legendre transform do not commute, it is in general not viable to compute 
corrections to the Hamiltonian and then Legendre transform to find an effective Lagrangian 
to read off possible Lorentz violations. The corrected Hamiltonian itself would not be of 
higher order in time derivatives and so the corresponding Lagrangian would be analytic 
in the expansion parameters. But if higher derivatives for the full expressions have to be 
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expected, the Legendre transform of the perturbed Hamiltonian would not coincide with 
the perturbed full Lagrangian. In particular if there are higher spatial derivatives in the 
Hamiltonian, as in any spatially discrete theory, there are two possibilities much more com- 
plicated to distinguish perturbatively: Either there are no higher time derivatives in the 
corresponding full Lagrangian, which would break Lorentz symmetry since there are higher 
space derivatives; or there are higher time derivatives, in which case the theory may or may 
not be Lorentz invariant. 

We illustrate these issues with an example for a discrete theory with finitely many degrees 
of freedom. Let the action be S = ^2 n eL(q n ,q n+ i) with Lagrangian L(q n ,q n+ i) = (q n+ i — 
q n ) 2 /2e 2 , a discretization of a free particle with discrete time step e. If we define the momen- 
tum by p n := edL/dq n+ i and the Hamiltonian by L(q n , q n+1 ) = p n {q n+ i - q n )/e- H(q n ,p n ), 
we obtain p n = (q n+1 - q n )/e and H(q n ,p n ) = p 2 n /2. 

We now assume that e is small and approximate the discrete values q n by a continuous 
function q(t) such that q n = q(en). Thus, 

2 



00 k-1 




' I J- 

2 



-(q' + eqq + e^qq^ + lq^ + O^)) 



,fe=i 

which yields a higher derivative theory at second order in e, which we will use in what 
follows. (To linear order in e, however, the theory is not higher derivative since the only 
correction is a total derivative.) Only q itself and the first derivative q are independent 
variables since q^ can be removed from the Lagrangian by integrating by parts. Removing 
all the total derivatives and higher orders in e results in the Lagrangian 

which after performing a (higher derivative) Legendre transform gives momenta 

dL d dL . 1 2 ,„s 
^ '= ^-dt^ = q+ ~2 e q 
dL 1 2 .. 

:= m = ~~i e q 



and the Hamiltonian 



H(q, 7i q , q, Tig) = qix q + qTTg — L = qn q — ^q 2 - 6e 2 n 2 . (16) 



As anticipated, the Hamiltonian is not analytic in e. 
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If we had started in a Hamiltonian formulation and gotten our discrete formulation there, 
as in loop quantum gravity, we would have proceeded differently. First, expanding a Hamil- 
tonian formulation does not introduce new degrees of freedom such as q above, which is 
independent of q and has its own momentum. Moreover, the symplectic structure would be 
left untouched since it is independent of the Hamiltonian (while a Lagrangian formulation 
mixes symplectic structure and dynamics). Thus, we would still work with the unperturbed 
momentum n q = q for the only degree of freedom q. This fact can also be derived from 
the Hamiltonian (|16|) by solving the Hamiltonian equation of motion q = dH/dir q for 7c q (q) 
which is now assumed not to be an independent variable. We obtain q = q — (q — n q )dq/ d7i q , 
which, since now q is by assumption no longer independent of 7r q , immediately gives ir q = q. 

The second difference is that only perturbative corrections to the Hamiltonian could 
appear, and not the 1/e-term above. If we remove this term and use the unperturbed 
momentum, we obtain H = \it 2 which coincides with the continuous approximation of 
the full discrete Hamiltonian. Thus, to this order of perturbation theory, the Hamiltonian 
would not show any corrections unlike the Lagrangian. The reason is that in this example 
the perturbative corrections are all of higher derivative form, which cannot be seen in the 
Hamiltonian. 

The example can easily be extended in such a way that also the Hamiltonian re- 
ceives perturbative corrections. We now use two discrete coordinates and Lagrangian 
L = ((<7m,n+i ~ qm,n) 2 — (<Zm+i,n — <?m,n) 2 )/2e 2 where we interpret m as a discrete space 
coordinate and n as discrete time, as above. The Lagrangian is symmetric under the ex- 
change of m with n, which mimics a space-time symmetry. The momentum of q m , n now is 
Pm,n = {q m ,n+i ~ Qm,n)/^ and the Hamiltonian H = \{j> 2 m .n + (<?m+i,n - q m ,n) 2 /^)- 

Perturbing the Lagrangian leads to L — \{q 2 — (q 1 ) 2 + e(qq — q'q") + e 2 (^qq^ + \q 2 — 
\q'q'" — |(g") 2 ) + 0(e 3 )). The perturbed momenta are as before, and the Legendre transform 
of the perturbed Lagrangian with total derivatives removed as before yields 

H(q, 7i q , q, 7r q ) = qn q - l -q 2 + ^{q' f + ^ 2 (q") 2 - Qe^ 2 . 

With the prescription above, a perturbation of the Hamiltonian would have led to the 
analytic expression H = ^(^ q + {q') 2 +j2€- 2 {q") 2 ), which only shows the higher order correction 
(which is of higher derivative form only in space but not in time) , but not the non-analytic 
correction coming from the higher derivative nature. 
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These examples have important hints for the calculation of corrected dispersion relations 
and the issue of Lorentz covariance. Since only higher order corrections will be seen when a 
Hamiltonian is perturbed, Lorentz violations are bound to appear as a consequence of this 
way of doing the calculation. Space and time derivatives of the classical fields have to be 
related in the Lagrangian in a way dictated by the symmetry. If those terms are torn apart, 
because one computes the Lagrangian from a perturbed Hamiltonian which only sees higher 
space derivatives but not higher time derivatives in its corrections, Lorentz invariance will 
be violated. This kind of violation of Lorentz symmetry is not a consequence of the theory 
but of the way to perform perturbative calculations. 

We present one more example showing the role of higher derivatives in Lorentz invariant 
theories. We use the Lagrangian 

L = ~J (■?/>(□ + eU 2 )i) + m 2 V> 2 ) 

for a scalar of mass m. It leads to the field equation 

-(□ + eD 2 )^ = m 2 ip 

which is Lorentz invariant. The dispersion relation can be computed from the plane wave 
ansatz ip(x,t) = exp(i(Et — kx)) and takes the form 

eE A + (1 - 2ek 2 )E 2 + k 2 (ek 2 - 1) = 

such that 

E 2 = k 2 — ± Vl + 4em 2 /2e . 

If e <C m~ 2 , we obtain 

E 2 = k 2 - -(1 =F (1 + 2em 2 - 2e 2 m 4 )) 

with two non-analytic solutions, which have to be discarded in a perturbative situation, and 
the corrected relation 

E 2 = k 2 + m 2 - em 4 + 0(e 2 ) 

which is Lorentz invariant. 

The Hamiltonian situation of this example is as follows. We have momenta ir^ = ip — 
2eAip + etp and = —eip leading to the Hamiltonian 

1 • • 1 1 • 1 
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Again, this is non-analytic in e, but would lead to Lorentz invariant equations of motion. 
If, on the other hand, we had started with perturbing a Hamiltonian, we could only have 
seen the analytic part and would not have introduced new degrees of freedom and instead 
used 7r^ = ip (this again also follows from the equations of motion under the assumption of 
having no additional degrees of freedom: ip = dH/dir^ = —ipdip/dn^ + ip + ^^dip/ii^ implies 
immediately = ip if, as per our assumption, dip/dit^ ^ 0). Thus, we would have arrived 
at a Hamiltonian 

H =o^i~ + + e ( 27r V> A7r </> + V'AV) 
and perturbed equation of motion 

■0 = Aip - m 2 ip + e(A 2 ip - 2m 2 Aip) . 

The dispersion relation for this equation is 

E 2 = k 2 + m 2 -ek 2 (k 2 + 2m 2 ) 

which does break Lorentz invariance. 



IV. CONCLUSIONS 

The observations presented here beg the question of what is the correct procedure to com- 
pute modified dispersion relations from a Hamiltonian point of view when higher derivative 
terms have to be expected. Such a procedure has to be amended incorporating the semi- 
classical dynamics in a more controlled way. This in particular has to take care of new 
degrees of freedom that emerge from higher derivative theories. One possibility is to derive 
a full, non-perturbative discrete Hamiltonian from the quantum theory, which is understood 
as a classical object but on a discrete space e.g. [3]. Before one expands and computes 
equations of motion, one has to transform to a Lagrangian, also on discrete space and time. 
From then on one can work with perturbative expansions and compute modified dispersion 
relations. 

There are obvious difficulties in the way of implementing this procedure in loop quantum 
gravity since already calculations with a perturbed Hamiltonian are cumbersome. It should 
however be kept in mind that the calculations done up to now (including the model of 
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the previous section) can only yield preliminary results and that a definite answer to the 
question of Lorentz violation by loop quantum gravity definitely has to await a more complete 
treatment, possibly along the lines sketched above. Alternatively, perturbative Hamiltonian 
techniques for effective actions, which also allow to see additional degrees of freedom coming 
from higher derivatives, can be developed. This approach, which is now under investigation, 
would allow to perform the perturbation expansion at the Hamiltonian level all the time. 

We expect that the model presented in the first part of this paper can be used for a 
first step in applying those methods to the issue of dispersion relations. As we showed, it 
shares most qualitative and even some quantitative features with more elaborate calculations 
and thus is simple but reliable. It can therefore play a role in deriving modified dispersion 
relations that better take into account the higher derivative nature. 
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